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We develop an explicit and tractable representation of a twist-grain-boundary phase of a smectic- 
A liquid crystal. This allows us to calculate the interaction energy between grain boundaries and 
the relative contributions from the bending and compression deformations. We discuss the special 
stability of the tt/2 grain boundaries and discuss the relation of this structure to the Schwarz D 
surface. 
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Topological defects are often the essential degrees of 
freedom. They are the focus in the study of some phase 
transitions Q, , high-temperature superconductors Q , 
and liquid crystalline phases Q. In the latter, there is 
a necessary connection between the topology of the de- 
fects and their geometry which is on the one hand the- 
oretically challenging while, on the other hand, exper- 
imentally accessible through real-space, freeze-fracture 
imaging The twist grain boundary (TGB) phase 

of smectic-A liquid crystals |6(, has an arrangement of 
screw dislocations which alter the geometry of the uni- 
form, flat layers into a discretely rotating layered struc- 
ture. Though topology constrains the geometry, it does 
not specify it. Rather, the free energy of the deformed 
smectic layers sets the periodicity of the lattice. Prior 
analysis |7| relied on linear elasticity to study small an- 
gle grain boundaries. However, it has been shown that 
when the angles and deformations are large, the energet- 
ics of the rotationally-invariant nonlinear theory are not 
only quantitatively, but qualitatively different than in the 
linear theory 0, 0, ■ In order to reconcile our under- 
standing of the linear theory with the nonlinear elastic- 
ity and in light of the recently observed large angle 
grain boundaries, here we develop a full nonlinear the- 
ory of the largest angle grain boundaries allowed, with 
rotations of 7r/2. To do this we explicitly sum the topo- 
logical defects to render a closed-form expression which 
is an exact solution of the linear elasticity theory. Un- 
like parametric representations of surfaces, our surface is 
given as a multi-valued height function which allows us 
to directly calculate the compression energy and allows 
tractable comparison with real space images. We directly 
calculate the energetics of space-filling TGB structures 
analytically and find that grain boundaries interact ex- 
ponentially with separation. 



Smectic order is characterized by a periodic mass den- 
sity p = po + pi cos [27r$(x)/a], where po and p\ are 
constant amplitudes, set at the nematic-to-smectic phase 
transition, and a is the equilibrium layer spacing. The 
smectic layers are defined via the level sets of $ through 
$(x)/a = n 6 Z, and the elastic free energy has two 
terms, 




FIG. 1: Schnerk's first surface, with charge +2 and —2 screw 
dislocations. Note that the —2 dislocations lie at the center 
of the rectangle made by the adjacent +2 dislocations. We 
choose 6 = tp and k 2 « -0.03033 so that K'(fc)/K(fc) = 2 - i. 
We are reminded that ". . . imaginary things are often easier 
to see than real ones." fl2]| 
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where B is the compression modulus, A 2 = Ki/B, K\ is 
the bending modulus, and H = • [V<&/|V$|] is the 
mean curvature of the layers. To study the energetics 
of defect-laden structures, we find an expression for the 
phase function <!> of the surface and use that to evaluate 
(|T|). For instance, a single screw dislocation at the origin 
is described by a helicoid |J| (here and throughout, x,y, 
and z are coordinates in R 3 and w = x + iy): 

Escrow = 7 2 - 7T~ arctan (-) = 72 - -^-Imlnw (2) 
27T \x/ 2ir 

where the Burgers scalar b/a must be an integer so that 
the mass density is single valued. The prefactor 7 is 
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necessary in order to allow |V$| — > 1 away from the 
defect. For the screw dislocation this results in 7 2 = 
1 — lim u ,^ 00 [&/(27r|w|)] 2 = 1. The helicoid is an extremal 
of 0} as well as an extremal of the often used quadratic 
free energy (which controls linear elasticity) 
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written in terms of u = z — <£> and Vj_ = xd x + yd y . 
Unlike the nonlinear theory, however, a screw dislocation 
has vanishing energy in the linear theory because $ scrcw 
is a harmonic function of w, i.e. V^^c-ow = 0. 

Do multiple screw dislocations interact? In the lin- 
ear theory J2J there is no interaction. The inclusion of 
non-Goldstone, nematic director modes results in long- 
distance exponential interactions as arise in flux lines [6| . 
However, here we consider the interactions arising from 
the necessary nonlinearities of a rotationally-invarient 
free energy Our study proceeds by considering so- 
lutions of the linear theory and calculating their energy 
in the nonlinear theory. Though these solutions are not 
minimizers of Q they have the desired topology of the 
smectic configurations of interest and, importantly, be- 
come exact in the limit of infinite separation between 
dislocations. We begin with a single twist-grain bound- 
ary, a row of dislocations along the cc-axis with uniform 
spacing £4 ■ The phase field is a sum of individual heli- 
coids: 
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Utilizing the infinite product sin u; = wj} n _^ (l 
follows that ^3 ( U P to unimportant constants) 
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Since the grain boundary lies on the y = axis, set- 
ting the compression strain to zero at y — ±00 com- 
pels us to set 7 2 = 1 — b 2 /(2£d) 2 . The layer normal of 
a single grain boundary, N = V$/|V$|, rotates from 
N_ = [-227)0,7] to N + = [217,0,7] as y goes from -00 
to 00, implying a uniform rotation of the layers across 
the grain boundary by an angle sin a = b^/id EIIH- 

Because a rotation by 7r amounts to no rotation at 
all, there is necessarily a dual description of a single 
grain boundary as a rotation of it — a by viewing the 
dislocations as parallel to the £-axis. This dual de- 
scription reflects the geometric nature of the defects 
and can be seen directly from the parametric equation 
for the level sets of $ r ow, tan(27T7z/6) tan(7rx/£d) = 
t&nh(iry /£d) 0] Thus, after a rotation around the j/-axis 
by 7r/2, (x,z) — > (z, — x), and the layer normal rotates 
from N_ = [7,0,^] to N+ = [-7,0,^7]. The rota- 
tion angle becomes sin a = — 67/^d, for a total rotation 



of 7r — a. Note that the sense of the rotation has been 
reversed or, equivalently, 6 — s- — 6 |Q] . 

This duality is the lynchpin for our further analysis. 
When piecing together single grain boundaries into a 
TGB phase, the defects in the grain boundaries must 
rotate with the smectic layers, pulling the topological 
defects along with the very geometry they create. In the 
case of 7r/2 grain boundaries a special simplification oc- 
curs. If the first grain boundary has defects along the 
z-axis, then the two adjacent boundaries should have 
defects along the x-axis. Employing the duality which 
swaps x and z, we can choose to view the adjacent bound- 
aries as made of defects along the z-axis with the oppo- 
site Burgers scalar. Thus, in the case of ir/2 boundaries, 
we have a structure with only parallel screw dislocations 
along the the z direction, with alternating signs! The 
phase field for the sum of individual grain boundaries is 
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where r is a complex number which generates the ap- 
propriate translations along the y direction. Throughout 
we set r = 2i£b/£d + 1, where £b is the distance between 
the grain boundaries p. Note that this sum alternates, 
reflecting the alternating sign of the defects in adjacent 
boundaries. 

The infinite sum in equation © can be put into closed 
form by observing that the exponential of the sum is 
doubly-periodic and, through rescaling of x and y, shares 
all the poles and zeroes of the Jacobi elliptic function 
or equivalently, through one of the estab- 
lished infinite products for sn(u, k). This sum will gen- 
erate a surface of the desired topology which is also a 
harmonic function and thus a minimizer of the quadratic 
free energy. We arrive at the exact summation of screw 
dislocations for a ir/2 TGB structure: 
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where 8 = 2K(k)/£d, ip = ReK'(fc)/£& are scale factors, 
K(fe) is the complete elliptic integral of the first kind, 
K'(/c) = K(-\/l — k 2 ), and k is the elliptic modulus, which 
for our purposes must be pure imaginary (i.e. k 2 is real 
and nonpositive) The level sets of ^tgb, dubbed 

Schnerk's first surface, are shown in Fig. ^fo r @b = £d- 
To simplify our notation, we define £ = 8x + itpy, use 
Glaisher's notation (cs for cn/sn, etc.) for the elliptic 
functions [l^. and suppress the elliptic modulus k. 

The orthorhombic symmetry of Schnerk's surface sug- 
gests that, like the single grain boundary, it can be viewed 
as one of three orthogonal arrangements of screw defects. 
To see this, we note that the level sets satisfy 
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We recognize sc(C)/dn(£) as the elliptic generalization of 
tan £. By considering the zeroes and poles of the elliptic 
functions, it is possible to view the surface as being com- 
posed of oppositely charged, staggered defects along x or 
y, instead of z, generalizing the duality of a single grain 
boundary 0- 

The closed-form expression for $tgb allows us to study 
the energetics of this structure by use of the established 
properties of elliptic functions. For instance, the com- 
pression strain u zz = [1 — (V$) 2 ]/2 is the somewhat 
cumbersome 



1 

u zz = - 



l- 7 2 --^(# 2 +^ 2 )|csCdnCr 



(9) 



This expression demonstrates that we can choose the pair 
(9, ip) in order to simplify our analysis; because we can 
set the periodicity of the lattice by altering (9, ip) or k, 
we can freely set the pair at our convenience. Though 
these are distinct deformations of the layered structure 
they share the same periodicity and topology. We will 
focus on the case ip = 9 in the following, though our re- 
sults do not change qualitatively for other choices. The 
elliptic modulus is set by ReK'(fc)/K(fc) = 2£ b /£ d . Re- 
call that in the case of a single grain boundary, we set 
7 by considering y — ±oo. Here, the structure is triply 
periodic and there is no "infinity" . We choose instead to 
have the compression vanish halfway between the grain 
boundaries, e.g. along y = £b/2 or x = £<j/4. Because 
these lines are where we would measure the rotation of 
the layers, this is a natural choice and, for k 2 < 0, the 
compression strain is constant along these lines. This al- 
lows us to set both 7 and l d - Though in principle we 
should choose 7 to minimize the compression energy for 
a single periodic domain, as ib/^d — * 00 the latter proce- 
dure will yield our solution. 

We make use of the following expansion [l6| in terms 
of q = exp [-ttK'/K] = - exp [-2ir£ b /£ d ] : 
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The utility of this expression is that it and its derivative 
have, as their leading terms, the trigonometric functions 
present in the single grain boundary (yEJ), allowing us to 
isolate the energetic corrections arising from interactions. 
To calculate the interaction energy between grain bound- 
aries, we expand the free energy in powers of q. Since 
q = — exp [—27r£b/£d], the interactions will fall off expo- 
nentially with £b, our central result. We note that this 
conclusion is independent of our choices of 9 and ip as 
long as we are in a regime where q is small. For large 
£b/£d, q is small making this a good expansion. Note 
that even for the symmetric case where £ b = id, shown 
in Fig. 1, we have q — —e~ 27r w —0.002, small enough 



to use l)10|l reliably. We compare this to the interaction 
between defects in a linear theory with a director field: 
in the latter case the decay length is A, not £ d . The inter- 
action in the nonlinear theory arises from long-distance 
strain not from the decay of gauge- like director modes. 

In order to find the interaction between grain bound- 
aries, we expand the energy per unit area in terms of 
q and remove the energy of L y /£b isolated grain bound- 
aries, where L y is the y dimension of the system. The 
compression interaction per area A is proportional to 
qj \nq and we find: 
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where C is a positive constant of order unity, L z is the 
^-dimension of the system, and an elastic cutoff length 
£ is introduced to cutoff a |w| divergence in u 2 zz near 
the origin. This cutoff is necessary in the case of a sin- 
gle dislocation and a single grain boundary || as well. 
Since q < 0, we find an attractive, exponential inter- 
action. The attraction comes as no surprise since the 
adjacent grain boundaries are made of opposite signed 
dislocations. Note that the constant term in gets a 
contribution from the "tails" of the single grain bound- 
aries and is independent of both the Burgers scalar and 
the cutoff. The numerical details of (fTTf) depend on the 
geometry of the cores, but our result is representative 
of the energetics. Minimizing over £ gives £ oc b as in 
nonlinear theories of edge dislocations 0, Q . 

The bending energy does not require a cutoff, as the 
mean curvature is finite everywhere (see Fig. |2J) . This en- 
ergy may also be expanded in powers of q and we would, 
again, find an exponential interaction. In the two special 
cases of ip = \/29 for any k and 9 — ip with k 2 = — 1, 
the interaction between grain boundaries is purely repul- 
sive, the first case recapitulating the vanishing of H for 
Scherk's first surface |8| . Whether the curvature interac- 
tion is always repulsive is under investigation |l7|. 

Our construction demonstrates that the it/ '2 TGB 
structure is, in fact, charge neutral from the point of 
view of the screw dislocations. Indeed, any TGB struc- 
ture with rotation angle a — n/n for n S Z will be charge 
neutral - the screw dislocations in one boundary will have 
equal and oppositely charged defects in the n th further 
grain boundary. However, the n/2 case is special because 
the cancellation occurs in the adjacent grain boundary 
which suggests that the 7r/2 structure is especially sta- 
ble. Though the structure we have studied here is not 
chiral, a nematic director can rotate uniformly through 
the structure, pulling away from the surface normal as 
necessary. For n = [cos(|j^), 0, sin(|^)], we find the av- 
erage alignment of the director with the layer normal to 
be S = §((N • n) 2 ) — 5 ~ 0.7 for the geometry shown in 
Fig. 1 . This value suggests that the director and the layer 
normal are more or less aligned. If the geometry of the 
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FIG. 2: The unit cell of Schnerk's surface shares the topology 
of the Schwarz D surface. However, Schnerk's surface is not 
minimal. We project the value of \H\, the magnitude of the 
mean curvature of the surface, in grayscale onto the floor 
(black is zero curvature, white is the maximum curvature). 



mesogens were to favor positive saddle-splay or, equiva- 
lently, negative Gaussian curvature then this would fur- 
ther stabilize this phase over flat layers. It may be pos- 
sible to study grain boundaries with rotation angles ir/n 
by grouping blocks together into "superblocks" which in- 
teract as effective ir/2 grain boundaries. Whether this 
can be made precise, even perturbatively in q, is under 
investigation |l7j ]. 

In closing, we note that the unit cell of Schnerk's first 
surface shown in Fig. [21 has the same bicontinuous topol- 
ogy as the minimal Schwarz D surface. In Fig. |2 we 
have also indicated the magnitude of the mean curva- 
ture, \H\, which does not vanish everywhere. We recall 
that a single grain boundary can be made minimal by a 
stretch of sec(a/2) along the y-axis to become Scherk's 
first surface [8|. A possible generalization, in the spirit of 
recent solutions of maximal surfaces in Minkowski space 
fl8| . is to replace tan(27T7z/6) with its elliptic analog 
sc(2ir'yz/b)/dn(2Trjz/b) in JSJ to achieve greater symme- 
try between the x,y, and z coordinates. Similarly, the 
parametric Weierstrafi-Enneper representation 0] of the 
coordinates of a minimal surface can, in some cases, be 
reduced to elliptic functions [2(J. We are unaware of 
choices of 6, tp, or k or other simple deformations that 
make the curvature of Schnerk's surface constant. 



We have directly summed the phase fields for an infi- 
nite array of screw dislocations which generates the topol- 
ogy of 7r/2 TGB phase and which is a solution of the 
linear elasticity theory. This provided an analytically 
tractable representation of the structure which we used 
to study the interaction energy between grain boundaries 
in a rotationally-invariant elasticity. We find exponential 
interactions - the true minimizer should have interactions 
at least as weak. Further work will explore additional 
deformations, the character of the cutoff £, and make 
intimate comparison with experiment . 
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